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Abstract The scaling of entanglement entropy for the nearest neighbor an- 
tiferromagnetic Heisenberg spin model is studied computationally for clusters 
joined by a single bond. Bisecting the balanced three legged Bethe Cluster, 
gives a second Renyi entropy and the valence bond entropy which scales as 
the number of sites in the cluster. For the analogous situation with square 
clusters, i.e. two L x L clusters joined by a single bond, numerical results 
suggest that the second Renyi entropy and the valence bond entropy scales 
as L. For both systems, the environment and the system are connected by 
the single bond and interaction is short range. The entropy is not constant 
with system size as suggested by the area law. 

Keywords Entanglement entropy ■ Area law • Valence Bond Monte Carlo 


1 Introduction 

This paper is a numerical investigation of entanglement entropy of the nearest 
neighbor isotropic Heisenberg spin 1/2 model on clusters, in particular, a 
Bethe cluster and two L x L square clusters joined by a single bond. We are 
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Fig. 1 14 site three branch Bethe cluster. 


studying the ground state quantum mechanical properties of these models 
and the coupling J is taken to be the same on every bond (including the 
single bond joining the clusters), that is the interaction is antiferromagnetic. 
Recall that for a large Bethe cluster and for the square lattice there is very 
good numerical evidence that there is long range antiferromagnetic order 
(suitably defined for a Bethe cluster), though to the best of our knowledge 
there is no proof. The numerical methods used are spin wave theory J], 
direct diagonalization and valence bond Monte Carlo mm- There is no sign 
problem associated with the models so Monte Carlo is an effective numerical 
method. The quantities calculated by valence bond Monte Carlo are the 
valence bond entropy [4lf5] and the n = 2 Renyi generalized entropy S 2 [S'- 
Both these quantities are straightforward (given the techniques in ref. [23 ) 
to calculate by valence bond Monte Carlo. Note that Si, the Von Neumann 
is not so easy to calculate, however, it is believed that the same essential 
physics is contained in S 2 (but see ref. j6j). Generically, we will refer to all 
these entropies as entanglement entropy. We only use ’’balanced” clusters, 
where the number of even sites is equal to the number of odd sites, thus the 
ground state has spin 0 [7] and the ground state can be represented by a 
superposition of valence bond states [8ll2] , 

In particular, for the present study one of the clusters we shall consider is 
the three branched bond centered Bethe clusters, see figure 1 for an illustra¬ 
tion of the 14 site 3 branched bond centered cluster [S] . As discussed in [5] , if 
you bisect such a cluster, a simple argument shows, the valence bond entropy 
must scale as the number of sites N. In contrast, as we will later numeri¬ 
cally demonstrate, the spin-wave technique gives an entropy proportional to 
logiV. Because of this apparent contradiction it is important to use an un¬ 
biased numerical method, valence bond Monte Carlo, to calculate the Renyi 
generalized entropy S 2 . Due to recent advances in numerical techniques j3llTTj 
such calculations can be done accurately for large clusters. A priori, however, 
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one does not know if the clusters one can treat are large enough to see the 
asymptotic scaling law. One of the objects of the current investigation is to 
determine if one can realize the asymptotic regime with existing numerical 
techniques and hardware. 

Why are we interested in the scaling law for entanglement entropy? 
Firstly, if the entanglement entropy scales as In N, DMRG (Density Matrix 
Renormalization Group) is an effective, unbiased method to calculate the 
ground state properties. That is, one expects the number of states needed 
to describe the DMRG blocks to go as e Sl ~ O(N) not say e N . A still out¬ 
standing issue in condensed matter physics is how Neel order is destroyed by 
the addition of holes and ultimately is transformed into superconductivity. 
Since the Bethe cluster has Neel order at 1/2 filling [T21[T5] , the Hubbard or 
t-J model on the Bethe cluster would be an effective way to study this issue 
assuming the DMRG blocks scale with the number of sites, not the number 
of states, in the blocks. 

Secondly, unlike in one dimension, the status of the area law for the 
entanglement entropy is not as clear M- What are the conditions on the 
Hamiltonian and the cluster for the validity of the area law for a non one 
dimensional system? Naively, since a single bond connects the two halves of 
the Bethe cluster, one would expect either a constant or a logarithmic depen¬ 
dence of entanglement entropy on system size. However, for non interacting 
fermions, one sees N dependent entropy for the Bethe cluster and y/N or 
log N (depending on boundary conditions) for square clusters separated by a 
single bond [Hi. Is this a pathological feature of non interacting systems? 

There have been a number of very interesting papers |151I16[[T71I18| where 
examples of models exhibiting large entanglement entropies or violating the 
area law, are investigated. To best of our knowledge, the work developed 
in these papers does not apply directly to the particular systems we con¬ 
sider. Broadly speaking, it seems in the above papers, the lattice or cluster 
is straightforward to experimentally realize, while the Hamiltonian is diffi¬ 
cult to realize in a practical situation. For our models, the Hamiltonian is 
physically realistic, however, a very large Bethe cluster is hard to realize in 
an experiment m- However, there is no such difficulty in realizing the two 
L x L cluster system. 


2 Bethe-Cluster 

2.1 Spin-Wave Approach 

To calculate the Von Neumann entropy Si and the second generalized en¬ 
tropy S '.2 from a spin-wave approximation we use the approach of [I: which 
is easily applied to the Bethe Cluster. Let us consider initially a different 
situation from that considered in DMRG, namely, we take the subsystem in¬ 
side the system. As an example, consider in Figure 1, a three site subsystem, 
consisting of the sites labelled 1,2,3 inside the 14 site bond centered cluster. 
Figure 2 is the Von Neumann entropy for subsystems of size 3 to 63 for the 
254 and 510 site clusters. 
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sites in interior subsystem 

Fig. 2 Si vs. sites in the interior subsystem. The red circles are for the 254 site 
cluster and the blue diamonds are for 510 sites. The curve is a linear fit to the 
points for the 510 site cluster. 


We see the Von Neumann entropy scales with the number of sites in 
the interior cluster assuming one is sufficiently far from the boundary of 
the cluster. Note for the 254 site cluster the 63 site interior cluster is quite 
far from the linear fit in figure 2 while for the 510 site cluster the 63 site 
interior cluster is on the linear fit. The linear scaling is consistent with what 
one expect from the area law, since for an interior cluster, the number of 
boundary points scales with the number of sites in the cluster. 

Let us now examine a situation of greater similarity to that encountered 
in the blocking procedure in DMRG. Take a bond centered cluster (figure 1) 
and pick the subsystem to be the left half of the cluster. 

This is done in figure 3, the Von Neumann entropy and S 2 are plotted 
vs. the logarithm of the system size. We see both quantities scale as the 
logarithm of the system size, similar to a one dimensional system. However, 
as previously mentioned, by the argument of ref. |5] , the valence bond entropy 
scales as the number of sites. Naively, one would expect, say the valence bond 
entropy and Si (or S 2 ) to scale the same way with system size ( at least up to 
logarithms EDI). Thus either the valence bond entropy scales differently from 
other entropies or the spin wave calculation gives an incorrect result. In figure 
3, the green squares refer to exact diagonalization results for S 2 for system 
sizes 6, 14, and 30. The system sizes accessible to direct diagonalization (the 
30 site cluster has state space of dimension approximately 1.6 x 10 8 ) are too 
small to infer the scaling of S 2 with system size. Nonetheless, some insight 
can be gained from plotting, in figure 4, the entanglement entropy for the 
30 site cluster vs. the number of sites in the interior subsystem. Perhaps 
not surprisingly, the entanglement entropy increases as the number of bonds 
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Fig. 3 Entanglement Entropy vs. ln(cluster size) for subsystems bisecting the sys¬ 
tem. The red circles are for Si while the blue diamonds are for S 2 . The green 
squares are exact diagonalization results for S 2 . 


connecting the subsystem to the system increases (up to 7 sites), after which 
the number of connections decreases and the entropy decreases. 


2.2 Valence Bond Monte Carlo 

We thus turn to valence bond Monte Carlo J2] as a method to compute, 
essentially exactly, the properties of the Heisenberg model on a Bethe Cluster. 
In figure 5, the valence bond entropy is plotted vs. system size for a bisected 
system. The valence bond entropy m is a natural quantity to compute 
with valence bond Monte Carlo, as the basis consists of valence bonds and 
the valence bond entropy counts the number of valence bonds leaving the 
subsystem. From the figure, one sees that the valence bond entropy scales 
with system size. Of course, given the argument in [9], this is no surprise. 
Further insight can be obtained by looking at the value of the valence bond 
entropy. By the argument of [9j, at least 1/3 of the valence bonds (for 1/2 
the cluster, i.e. for a 1022 site cluster, roughly 170 bonds) must connect the 
two halves of the cluster; from figure 3, we see for large clusters, very close 
to 1/3 of the bonds connect the two halves. 

Let us now consider the n = 2 Renyi generalized entropy S- 2 - Due to 
recent advances in computational technique this quantity can be calculated 
with valence bond Monte Carlo. We apply the methods developed in ref. [2, 
see also DU; S%{pa) is calculated as — In ((Swap A )) where A is the subsystem 
and Swap A is a swap operator. (Swap A ) is calculated in the simplest formu- 
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sites in interior subsystem 

Fig. 4 Entanglement Entropy vs. number of sites in the interior subsystem for a 
30 site cluster. The red circles are for Si while the green squares are for S 2 . The 
points were calculated by exact diagonalization. 
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Fig. 5 Valence Bond Entropy vs. system size for a bisected system. 
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sites in subsystem 

Fig. 6 S 2 vs. sites in the subsystem for a 254 site cluster. The red crosses are 
calculated from the r = 1 ratio method, while the blue diamonds use the shell 
technique. Statistical errors are smaller than the symbols in the figure. 


lation by a double projection Monte Carlo algorithm. In a more sophisticated 
approach the ratios 


(Swap A i+r) 

(Swap A i) 

are computed from Monte Carlo; from these ratios (Swap A ) is then cal¬ 
culated. Here r + i is a symbolic notation for a subsystem bigger than the 
subsystem i. For system sizes 6,14, 30 and 62 sites we have used three ’’differ¬ 
ent” approaches: the naive (no ratio) approach, the sophisticated approach 
where i + r consists of the next shell in the Bethe cluster and the brute force 
’’sophisticated” approach where ”r” is only one site. Recall there is a shell 
or layer structure for Bethe clusters, i.e. for the 30 site cluster (take 1/2 the 
cluster) the first layer has 1 site, layer 2 has 2 sites, layer 3 has 4 sites, layer 4 
has 8 sites. All three approaches agree, to within the statistical errors. Exact 
diagonalization results for 6,14 and 30 site clusters are also in agreement with 
the Monte Carlo calculations. 

In figure 6, we plot S 2 vs sites in the subsystem for a 254 site cluster. The 
subsystem is taken to be an interior subsystem as in figure 4. The red crosses 
are calculated from the r=l ratio method, while the blue diamonds use the 
shell technique. Statistical errors are smaller than the symbols in the figure. 
We see that S 2 appears to be a continuous piece wise linear function of the 
number of sites in the subsystem. It is linear within the shell with a kink in 
going from one shell to another. The kink (discontinuity in the derivative) is 
smaller in the shells near the center of the cluster. The decrease in entropy 
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sites in interior subsystem 

Fig. 7 S 2 vs. sites in the subsystem. Green squares are for 126 site system, red 
circles are for 254 site system and blue diamonds are for 510 site system. Statistical 
errors are smaller than the symbols in the figure. 


for the final shell is presumably caused by the spins on the boundary of the 
cluster that are no longer connected by bonds to spins outside the subsystem. 

We next consider a situation similar to figure 2. In figure 7, we take 
an interior subsystem and calculate S 2 via valence bond Monte Carlo for 
systems of size 126 sites , green squares, 254 sites, red circles and 510 sites, 
blue diamonds. Error bars would be smaller than the symbols in the figure. 
It appears that S -2 scales with number of sites in the interior subsystem for a 
sufficiently small subsystem relative to the system as expected from the area 
law. 

Finally in figure 8, we plot S 2 vs. system sizes for subsystems that consist 
of half the system. The blue diamonds are the spin wave results, while the 
red circles are calculated by valence bond Monte Carlo, r=l. The green cross 
is for the shell valence Monte Carlo method but breaking up the large shells 
to 32 sites. For system sizes less than or equal to 62, at least qualitatively, 
Monte Carlo and spin wave agrees; however as the system sizes grow past 
this point there is an increasing discrepancy. It appears the spin wave result 
increases logarithmically with the number of sites while the Monte Carlo 
results give linear dependence on the number of sites in the system. Does 
this crossover with system size make sense? 

One can present a simple minded argument that rationalizes this behav¬ 
ior. The relevant ’’aspect” ratio a is the ratio of the number of generations 
to the number of boundary points, for example for the 6 site cluster a = |, 
for 14 sites a = | for 30 sites a = | etc. . One would expect a must be small 
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sites in system 

Fig. 8 S 2 vs. sites in the system for subsystems bisecting the system. The blue 
diamonds are the spin wave results, while the red circles are calculated by valence 
bond Monte Carlo, r = 1. The green cross is for the shell method but breaking up 
the large shells to 32 sites. 


to see Bethe cluster rather than 1 dimensional behavior. Hence one would 
anticipate one needs to study system sizes substantially larger than 30 sites. 

A heuristic argument can be made for the magnitude of the slope of S 2 
vs. system size. Recall that S2 = — ln(0|Swap A |<^) where (j> = |V’o)|V , o) he. the 
replicated ground state. Note that |i/>o) = YYaf^ l a ) where |a) is a valence 
bond state and f a > 0. Then the normalization 

i = (m= E ( 2 ) 

oi,/3,oi f ,(3' 

and 

(^|Swap A |^>) = (3) 

E fafpfa’fp' (/3'| (a'ISwapA |a) |/3) (4) 

a,(3,a' ,/ 3 ' 

All the valence bond states have at least N/6 bonds connecting the two 
halves of the Bethe cluster. Typically, sites in A and an equal number 

in B are connected by valence bonds in both |a) and | /?). Assume states 
{P'\ and ( a'\ with the same valence bonds dominate the normalization and 
the expectation value of Swap A . The term (/3 , |(a'|Swap A |a)|/3) is reduced in 
comparison to {ft\{a'\a)\(3) by a factor of 2'Tn due to these bonds. Hence 

N N 

S 2 « -In2-B* = — In 2 
54 


( 5 ) 
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Fig. 9 4X4-4X4 cluster 


i.e. S 2 is proportional to N with a slope of about 1/78. This is comparison 
to the numerical results in figure 8 which give a slope of roughly yAj. 

3 Two Square Clusters Joined by a Single Bond 

As in @] , consider two Lx L clusters linked by a single bond. Here the bond 
is a spin-spin interaction of strength </, the same strength as all the other 
spin-spin interactions. The single bond is chosen to be as close to the middle 
of the side as possible. For L even (odd) this bond joins the L /2 ( (L + l)/2) 
site on the side facing the other cluster with the corresponding site on the 
other cluster (figure 9). There are several interesting aspects of this model. 
Firstly, this model can be realized in a physical system, either solid state or 
cold atom. Secondly, for L even, there exists valence bond patterns where 
no valence bond need to go from one L x L cluster to the other. Lastly, for 
non interacting fermions, at 1/2 filling, the entanglement entropy scales as 
L , not as a constant or In L. This is consistent with the Bethe cluster results, 
both for free fermions and the Heisenberg model, in that the ’’perimeter” of 
the model (the sites with only one bond) scales as the number of sites in the 
cluster. 

Figure 10 is a plot of the second Renyi entropy S 2 vs. L for L even taking 
the subsystem as half the cluster calculated using valence bond Monte Carlo. 
The values of S 2 are noticeably smaller than for comparably sized Bethe 
Clusters. However, S 2 does not appear to be a constant or depend on In L. 
The calculation appears, to within the large statistical errors and the limited 
system sizes, to be consistent again with S 2 scaling as L. To make the case 
for L more convincing, figure 11 is a plot of valence bond entropy vs. L for 
various system sizes and L even. One sees a rather clear indication of the 
valence bond entropy scaling with L. The same quantities are plotted for L 
odd in figure 12. Again one sees the valence bond entropy scaling as L. The 
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Fig. 10 S 2 vs L for L x L - L x L clusters, L even 
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Fig. 11 Valence Bond Entropy vs. system size for a bisected system, L even. Error 
bars are smaller than the data symbols. 



structure in the graph can be rationalized by noting that for L = 2n + 1 
n even, there is a valence bond pattern with only one valence bond joining 
the two L x L squares and that bond can be chosen to be the middle bond 
(where the interaction joining the squares is located). 


4 Conclusions 

Numerical evidence has been presented that for a Heisenberg model on clus¬ 
ters joined by a single interaction, the entanglement entropy scales as the 
perimeter, not as a constant, as suggested by the area law. This is consis¬ 
tent with work on non interacting fermions and demonstrates these earlier 
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Fig. 12 Valence Bond Entropy vs. system size for a bisected system, L odd. 

results are not an artifact of non interacting particles. As discussed in the in¬ 
troduction, this has ramifications for DMRG calculations on Bethe clusters. 
More generally, the scaling with perimeter places fundamental limitations on 
DMRG algorithms involving partitioning into clusters. This is analogous to 
the restrictions placed on heat engines by the thermodynamic entropy where 
” heat engine” in this case is a quantum calculation (done on a classical com¬ 
puter) and ’’entropy” is the entanglement entropy. In addition, two L x L 
clusters joined by a link is an experimentally realizable system. To directly 
measure entanglement entropy is difficult (though not impossible [21] ). How¬ 
ever, recent work on the area law suggests the obstruction to the area law 
in two dimensions ( see [22] p34-39 ) is the existence of edge states. Such 
(hypothetical) edge states may be easier to measure (and analyze) than di¬ 
rectly measuring the entanglement entropy. As well as telling us what we 
cannot do, insights from the second law of thermodynamics can be used for 
more rational engine design. An analogous situation could conceivably hold 
for the ’’area” law, where now the device analogous to a heat engine is a 
quantum computer. Proposed modular quantum computer architectures[231 
[M! are quite similar to the clusters joined by links considered in the present 
work. 
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